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O ■ Abstract 

7'- 

^ • In this work are computed analytical solutions for orbital motion on a background 

i_^2j. described by an Expanding Locally Anisotropic (ELA) metric ansatz. This met- 

ric interpolates between the Schwarzschild metric near the central mass and the 
Robertson- Walker metric describing the expanding cosmological background far 

T^ . from the central mass allowing for a fine-tuneable covariant parameterization of 

''^ I gravitational interactions corrections in between these two asymptotic limits. As- 

suming a non-varying gravitational constant, G = 0, it is discussed the variation 
of the Astronomical Unit (AU) obtained from numerical analysis of the Solar Sys- 

f-^ I tem dynamics, being shown that the corrections to the orbital periods on the So- 

CN [ lar System due to the decrease of the Sun's mass by radiation emission plus the 

General Relativity corrections due to the ELA metric background with respect to 
Schwarzschild backgrounds can be mapped to the reported yearly increase of the AU 

^ I through the corrections to Kepler's third law. Based on the value of the heuristic 

H ■ fit to the parameter AU corresponding to the more recent ephemerides of the Solar 

- - - System are derived bounds for the value of a constant parameter oq for the ELA 

metric as well as the maximal corrections to orbital precession and orbital radius 
variation within this framework. Hence it is shown that employing the ELA metric 
as a functional covariant parameterization to model gravitational interactions cor- 
rections within the Solar system allows to avoid the need for a varying AU and/or 
varying gravitational constant G. 



Keywords: orbital motion, Solar System, Astronomical Unit, local anisotropy, cosmological ex- 
pansion, variation of fundamental constants 



1 Introduction 

For numerical analysis of the Solar System dynamics it is chosen as standard of measurement 
the Astronomical Unit (AU). This standard of measurement fixes the measurement projection for 
spatial-lengths and, through Kepler's third law, relates the measurement projection for spatial- 
lengths with the measurement projection for temporal- lengths. Specifically the definition of the 
AU is stated in IS units as [IHl] 




AU = [GMq { -^^ = 149597870700m , (1) 



where Mq = 1.9891 x lO^'' kg is the Sun's mass, G = 6.67 x 10^^^ kg^^ m^ s^^ is the Gravitational 
constant, Tiau.o = 31562889.928... s is the Keplerian orbital period for a point mass in an 
elliptic orbit with semi-major axis of value 1 AU . This definition exactly matches Kepler's third 
law for a planet orbiting the Sun in an elliptic orbit with semi-major axis of length ri.orb = 1 ^U . 
This definition is based on the existence of a Keplerian constant of motion, specifically the angular 
momentum Jq (ri.orb = 1 ^U) = —y/GMAU (1 — e^), where e is the orbit eccentricity, hence time 
and space measurements are related by this classical conservation law stated in the definition ([1]). 

Considering experimental data from optical and range measurements of planetary orbital motion 
plus range measurement from orbiters and landers it was concluded that a yearly variation of 
the Astronomical Unit is required to fit the ephemerides in the Solar System. Hence adding the 
adicional parameter AU to the dynamical model, the best fit to this parameter was originally 
estimated to be [2l[3] 

(au) = +{15 ±4) cm yr-^ . (2) 

\ / Fit 

The value for this heuristic fit does not necessarily correspond to a physical radial distance vari- 
ation in the Solar System, instead it corresponds to a variation of the standard of measurement 
considered for numerical analysis of the Solar System ([T|) . 

As an alternative to the variation of the AU, it has been suggested that a variation of the 
product of the Gravitational constant by the Solar mass {GMq) would have a similar effect 
in the planetary dynamics. Specifically the best fit corresponds to the value {GMq)/{GMq) = 
— (5 zb 4) X 10~"'^^yr~^ |4J. We note that the variation of (GMq) has the effect of increasing 
both the orbital radius and the orbital period such that angular momentum is approximately 
conserved. Specifically the period varies as Tq/Tq « —2{GMq)/{GMq) and the orbital radius as 
AU/AU ^ -{G'Mq)/{GMq) [F,^. Hence, from Kepler's third law ([I]), a variation for (GMq) can 
be directly mapped to a non-null value for the parameter AU while keeping (GMq) [3] 

Au\ = 0.75 ±0.60 cm yr'^ . (3) 

The discrepancy between both fitted values ([2]) and ([3j) is due to the choice of the independent set 
of parameters considered in the fit. While in [2] both AU and AU are fitted, in [3] only (GMq) 
is fitted. This choice is justified in |4] by noting that AU and AU are correlated by 98.1%. We 
remark that it is also recognized in [[4j that the choice of whether to fit the parameter GM or AU 
is a matter of convenience and that both cannot be fitted simultaneously as the obtained value of 
AU does not exceeds the respective formal error. Actually both these choices are interpreted as 
a varying scaling of the measurement standard ([1]), hence a choice of the measurement projection 
(see [7] for further details), being directly mapped into eachother through Kepler's third law. 

This discussion is enough to conclude that some sort of unmodelled gravitational interaction seems 
to be acting in the Solar System which is effectively accounted for by an heuristic rescaling of 



the relation between spatial lengths and temporal lengths of the dynamical system model. The 
exact identification of the origin of such interaction is further difficulted due to the only available 
solution relying on the heuristic fit of one single isotropic parameter AU (or equivalently GM) 
across the Solar System. In this work our main objective is to employ a background described 
by an expanding locally anisotropic (EL A) metric [HI [9] to analytically model the corrections 
to orbital motion with respect to orbital motion on Schwarzschild background within the Solar 
System. In particular, based in analytical orbital solutions for the background described by the 
ELA metric we will explicitly compute the respective corrections to the orbital period due to two 
distinct contributions, the Sun's mass decrease corrections and the General Relativity corrections 
with respect to Schwarzschild backgrounds [10]. It is shown that it is possible to map these 
corrections to the heuristic variation of the AU, having simultaneously negligible contributions to 
the remaining orbital parameters. Considering the most reliable estimative for the average value 
of the parameter AU ([3]) we will compute bounds for the values of the metric functional parameter 
which, in turn, will allow to estimate the contributions due to the ELA metric background to the 
orbital period, orbital precession and orbital radius variation of the several planets within the 
Solar System. 

The ansatz for the ELA metric was originally suggested as a possible description of local matter 

distributions in the expanding universe |lipi2j. hence interpolating between the Schwarzschild 

(SC) metric [13] near the central mass and the Robertson- Walker (RW) metric [H] describing the 

expanding cosmological background far from the central mass. Hence the ELA metric generalizes 

the isotropic McVittie metric [15j and the anisotropic metrics considered in J16j , having the novelty 

of maintaining the SC event horizon and maintain as the only space-time singularity the SC mass 

pole at the origin such that the value of the Schwarzschild mass pole is maintained [HIT?]. This 

metric is locally anisotropic consistently with astrophysical observations [18'|19j converging at large 

radius to the RW metric such that global isotropy is maintained. Specifically the line-element for 

the ELA metric is 

ds'^ = {1 - Use) c'^ dt^ - rf {de^ + sin^ 9dip'^) 

(4) 



dri-/7— (1-C/sc) 2 + 2 cdt 



1 - Use ^ c 

where H = a/ a is the time dependent Hubble rate defined by the rate of change of the universe 
scale factor a. Use = 2GM/{c^ri) is the usual Schwarzschild gravitational potential, G is the 
Gravitational constant, M is the value of the Schwarzschild mass pole for the central mass being 
considered and c is the speed of light in vacuum. The exponent a is, generally, a function of the 
radial coordinate, here we will take the simplified ansatz studied in [8] 

2GM 
a(ri) = (ao - ai) + ai t/sc(n) = ("o - "i) + "i -^ — , (5) 

where oq and ai are numerical coefficients. The bound oq > 3 ensures that the SC radius 
'^i.sc = 2GM/(? is an event horizon and space-time is singularity free at this horizon, while for 
the bound ag > 5 space-time is asymptotically Ricci flat near the event horizon such that the 
SC metric is a good approximation in a neighbourhood of the point-like mass M. The coefficient 
«! < ensures that the singularity at the origin coincides with the SC mass-pole. In the following 
analysis we will consider it arbitrarily close to null, — 1 <C ai < 0, such that outside the event 
horizon its effects are negligible, hence a w qq for planetary orbits. For further details in the 
derivation of the ELA metric ansatz see [8] . 

We note that so far no direct physical interpretation for the functional parameter a exists, hence 
it is at most considered to be a functional parameterization of corrections to gravitational in- 
teractions for intermediate spatial-length scales. In particular we may assume that it is varying 



over the radial coordinate allowing to phenomenologically fit existing data hoping that the results 
obtained may shed some light on its physical interpretation (see for instance [2U]). In order to fit 
the constant oq to planetary motion within the Solar System maintaining the asymptotic limits 
of the ELA metric we are considering the simpler possible ansatz for this constant by considering 
that near the event horizon of the central mass (the Smi) the exponent is some constant ao.o ^ 3 
and that far from the event horizon (hence for all planetary orbits) the constant ao is approxi- 
mately given by a constant oq which can be either positive, either negative. Specifically we are 
considering the simplified ansatz 



ao.o > 3 , ri ~ Rqsc 



Qo 



(6) 



ao 



, n > Rqsc 



where Rqsc = 2GMq/c'^ ~ 2952.22 tti is the Sun's Schwarzschild radius and Mq is the Sun's 
mass. Based in this ansatz we will compute the analytical corrections to two body heliocentric 
orbital motion within the solar system on the ELA metric background (j3|) with respect to or- 
bital motion on Schwarzschild backgrounds |2HI22j. More accurate results including the many 
body gravitational interactions in the Solar System can only be computed by extensive numerical 
analysis including 10000 known bodies in the solar system |23j . These analysis are usually car- 
ried either in the Post-Newtonian formalism, either in the Post-Post-Newtonian formalism which 
includes both the General Relativity corrections to the classical Newton law of gravitation on 
Schwarzschild backgrounds, as well as corrections of extended theories of gravity such as Brans- 
Dicke gravity [24| . 

When required, for numerical evaluation of the Hubble rate H and the deceleration factor q 



of today's universe, we are considering the values Hq 

-2 



H\ 



t=to 



2.28 X 10^ 



-18 



and go 



a \a/ 



t=to 



0.582 |19j . As for the planetary orbital parameters considered for numerical 



evaluations we are considering the data presented in table [TJ 



Planet 


ri.orb (x10"to) 


e 


m{xW^kg) 


(Jorb (deg) 


<5pcr (deg) 


Mercury 


0.579091768 


0.20563069 


0.3302 


7.00487 


77.45645 


Venus 


1.08208926 


0.00677323 


4.8685 


3.39471 


131.53298 


Earth 


1.49597887 


0.01671022 


5.9736 





102.94719 


Mars 


2.27936637 


0.09341233 


0.64185 


1.85061 


336.04084 


Jupiter 


7.78412027 


0.04839266 


1898.6 


1.30530 


14.75385 


Saturn 


14.26725413 


0.05415060 


568.46 


2.48446 


92.43194 


Uranus 


28.70972220 


0.04716771 


86.832 


0.76986 


170.96424 


Neptune 


44.98252911 


0.00858587 


102.43 


1.76917 


44.97135 


Pluto 


59.06376272 


0.24880766 


0.0125 


17.14175 


224.06676 



Table 1: Planetary orbits parameters considered: the semi-major axis ri.orb, the eccentricity e, the mass, 
the orbital inclination Sorb o,nd the longitude of perihelion Sper 



This work is organized as follows. In section [2] are computed the analytical solutions for orbital 
motion on the ELA metric background being derived the General Relativity corrections to orbital 
precession and orbital period for such backgrounds. This derivation is carried considering a 
static elliptical orbit approximation. In section [3] are analysed circular orbits on the ELA metric 
background (jl]) and computed the orbital radius variation within this approximation. In section U] 
are analysed the corrections to Kepler's third law due to the decrease of the Sun's mass and due 
to the ELA metric background. In particular are derived estimative for the several contributions 
that allow to match the heuristic variation of the AU and computed the value for the ELA metric 



parameter that corresponds to the fitted value of the parameter AU ^, as well as the respective 
corrections to orbital precession and orbital radius variation. In the conclusions we shortly resume 
and discuss the results obtained in this work. 



2 Perturbative Static Elliptical Orbit Solutions 

In this section we derive analytical solutions for elliptical orbits on the background described 
by the ELA metric ansatz ^. In particular we will explicitly compute the General Relativity 
corrections on such backgrounds to the Keplerian elliptical orbit solution ri((p) = l/uQ^ip) with 

1 + e COs{ip) 2^ 1^-, 

uoW = -j , a = ri.orb(l - e ) , (7) 

where e is the orbit eccentricity and ri.orb is the orbital radius. We note that such corrections 
will generally include the General Relativity corrections on Schwarzschild backgrounds as well as 
corrections depending on the Hubble rate H. 

It is hard, if not impossible, to obtain an exact analytical solution considering the differential 
equations for a time varying Hubble rate H. The main difficulty is that energy conservation is 
no-longer given by a constant of motion, instead we have a non-linear second order differential 
equation on the function t coupled to the differential equation for ri. Hence, for technical simpli- 
fication purposes, we are taking the static orbit approximation by considering a fixed Hubble rate 
H = Hq corresponding to today's measured value for this rate. 

From the metric line-element (jl| let us consider the Lagrangian definition |21ll22j to order Hq 



C [^ ^^ ( ^^ ri\2,^ ^^ ^^\ / dt\ 






(8) 



-2i.o^(l-^sc)^-c^^-^3^^^J -^Hit ^^(^°^) 



This Lagrangian is a constant LIm = c and it is considered that the orbit of the test body 
is lying in the plane of constant coordinate 6 = tt/2 such that dO = and sin 9 = 1 |21il22j . 
The Lagrangian is independent of the coordinate (p, hence a constant of motion corresponding to 
angular momentum exists being given by the variational derivation of the Lagrangian with respect 
to dif/dr, 

J=— — = -r2^ (9) 

2mS^ ^ dr ' ^ ' 

dr 

Also, due to the Lagrangian ([8]) not depending explicitly on the time coordinate, a conserved 
constant of motion corresponding to energy exists being given by the functional variation of the 
Lagrangian with respect to c dt/dr 

2Eh 1 dC 



mc m 5{c dt/dr) 

This equation can be solved for c dt/dr such that replacing the obtained solution in the La- 
grangian ([5]), expressing the derivatives with respect to proper time dri/dr by the derivatives 



with respect to (/?, dri/dr = dri/dip x dip/dr and considering the change of variables u = 1/r, 
further differentiating with respect to y? and factoring out an overall factor of 2u'J^ (with the 
primed quantities representing derivation with respect to if), we obtain the approximate differen- 
tial equation of order Hq for the function u{if) describing an orbiting test mass in the gravitational 
field of a point-like central mass M 

GM 3GM 2 

GM /HoY f 2GM \-i+"o 

-"°^V~J V~^V (11) 

HoV 1 / 2GM \-i+°o / 2GM uqGM 

The terms in the first line match the usual terms obtained for Schwarzschild backgrounds and the 
terms in the second and third lines are the corrections due to the ELA metric background. 

We note that, although maintaining the terms of order Hq in the Lagrangian ([8]) which do not 
depend explicitly on the time coordinate we have neglected one term containing the factor Hq t. 
Explicitly it is the term —2qoHQtri (1 — Uscr°~^^'^{dt/dT) [dri/dr). Comparing the terms of 
order Hq in the Lagrangian with this term we conclude that this is a valid approximation as long 
as the value of the time coordinate is below the following bound 



1 ri ^2o_i Eh t„i3 

go dri/dr mc^ 



t<. —^{1-Usc)~~^—^^W^^ years. (12) 



This bound is well above any astrophysical measurement time span and has been obtained by 
considering the following simplified assumptions, within the solar system, from the experimental 
upper bounds on the orbital radius variations within the Solar System [7] we cosider the estimative 
for the ratio |ri/ri| > lO'^'', assume weak gravitational field Use ^ 1 fmd values of uq for which 
the approximation (1 — Uscr"°~^^''^ ~ 1 is valid. We note that this approximation will no longer 
be valid for very large values of the metric exponent, |qo| ^» 0. 

Noting that for orbits in the solar system the function u has relatively small values (0.5 x 10^^^ < 
u < 0.5 X 10~^'^ ?TT,~^, where ri orb is the orbit semi-major axis), with the objective of further 
simplifying the differential equation (jlip . we consider a series expansion on u of the terms of 
order Hq which is equivalent to an expansion on the weak gravitational field. Specifically, for a 
generic exponent p, the factor (1 — C/sc)^ has the following series expansion (l — 2GM u/c^) = 1 — 

p2GMu/c^+p{p-l) {2GM/c^f u'^/2-p{p-l){p-2) {2GM/c^f u^/Q+O (p"^ {2GM/c^f uA. 

We note that the full series is strictly convergent independently of the value of the exponent p 
as long as u < c^/(2GM), however an approximation to first order on u will only be valid as 
long as 2p GMjc? n < 1 , otherwise it is required to consider higher order terms to attain a valid 
approximation. Hence the differential equation (jlip is, to order u^, rewritten as 

n"(v.) + An(v.) « ^B + ^Cn^ 

HqY 1 a^Y GM 1 



J / u 



(H^y GM I ^.3. 



We note that with respect to the General Relativity orbit's equation on Schwarzschild back- 
grounds, there are the extra multiplicative factors A = 1 + 5a^B = \ + 5b and C = l + 6c- These 



factors differ from unity by the following additive constants 



Sc 



-«o (^ j (j^ + ("0 - 3ao + 2) 



2/9 ^ (GMHq 
■ o«o ("o - 3ao + 2j I 3 — 



2{GMf 
3c2 



1 + (ag - 7ao + 12) 



(14) 



{GMf 



So far we have not specify for which values of the coefficient uq the static approximation considered 
is valid. By comparing the leading order terms with the next to leading order terms we conclude 
that this perturbative equation is valid only for absolute values of the parameter uq up to 



|ao| < OO.max.pert 



2GM 



(15) 



Above this value it is either necessary to consider higher order terms of the series expansion or to 
consider the exact expressions. Nevertheless we remark that, for a fixed positive value of the radial 
coordinate ri, and larger positive values of the parameter uq > ao.max.pert the corrections given 
by the exact expression due to the ELA metric background will decrease significantly in absolute 
value becoming, for very large values of the parameter oq ^> ao.max.pert) negligible, while for larger 
negative values of the parameter ao < —ao.max.pert the corrections become more significant being 
unbounded from below. Hence, although for positive values of ao the approximation ()13p subject 
to the bound (llSp allows to establish a fairly good estimative for the maximum contribution of 
the corrections on the ELA metric background, for negative values of ao no bounds can be set 
for such contribution. In figure [Dare plotted the values of the exact and perturbative correction 
terms of order Hq in the differential equation (jlip . 
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Figure 1: Plot o] the exact (dashed line) and perturbative (continuous line) expressions for the corrections 
to the orbital differential equation ill]) on the ELA metric background as a function of the parameter ao 
for Earth's orbit. The perturbative regime is valid for |ao| < ao.max.pert ~ 5 x 10^ US]) : 

(a) plot of the exact expressions for aQ > 0, the corrections asymptotically vanish for large ao 3> ao.max.pert/ 

(b) plot of the exact and perturbative expressions for ao e] — 10^, 10'^[, the perturbative and exact expressions 
approximately match up to |ao| = ao.max.pert- 



In addition, with respect to the bound (fT2|) . we note that it is also obeyed as long as the bound (fT5|) 
is obeyed, hence for larger values of the coefficient ao the terms of order Hq explicitly depending 



on the time coordinate become relevant and must be included in the Lagrangian ^. For these 
cases there is no constant of motion directly associated with energy conservation. Specifically, 
conserved energy, would be given by the constant E = J dT{6C/6t — d/dT{6C/6t)), such that, due 
to the complexity of the full equations of motion, it would be preferable to consider a numerical 
analysis to compute orbital motion. 

We are proceeding assuming that the upper bound (jlSp is obeyed. To solve the differential 
equation (J13p we start by solving the differential equation considering only the dominant term in 
the right-hand side of p^ . hence obtaining [2HI22] 

i^o.m + ^'^o.H^ = ^p~ ^ =^ '^o.m = ^ , (16) 

where e is the elliptical orbit eccentricity and d is defined in terms of the ellipse semi-major axis 

'"l.orb as 

d = ri.orb(l-e') . (17) 

The standard General Relativity angular momentum Jq and the angular momentum J are ex- 
pressed in terms of the parameter d as 



Jo = -VGMd , 



B _ 1 
34 ^2 



(18) 



J = -^lGMd-f^-Joi5B-SA)j=j, , 



where to evaluate A and B, we have approximated the angular momentum by the respective 
Keplerian quantity, J ^ Jq. 

Next let us compute the corrections to the solution Uq jj2 by considering the remaining terms in 
the right-hand side of the differential equation (jl3|) evaluated for the function Uq jj2 (|16|) such that 
the full solution is 

U = Uqh^ + Uq^h^ + Uh^ . (19) 

Here the functions Ugr.h^ and Ujj2 correspond respectively to the corrections to the Keplerian 
orbit's solution due to the Schwarzschild background and due to the EL A metric background 
approximated to order Hq being, respectively, the solutions of the following differential equations 

// , _ 3GM 2 _3GM {l + ecos{VAip))^ 

2 



,, ^ , jHoV 1 ^ HoV GM 



J J Uqjj2 \ J J C Uqjj2 

HnV 1 



(20) 



+"0 ^- 



Jj {l + ecos{^/A^p))^ 



HoV GM 1 



J J c2 {l + ecos{^/Aip)Y 



such that we obtain 



Uqr.h^ = ^-^ ( (l + yj - — cos(2\/Ac^) + \/]4ev3sm(\/I(/9)j 
d^H^ faoGM {-A + e^) + 36"^ cos{2^/Aip) 



AJ2(l-e2)\^ c'^d 4(l-e2)(l + ecos(^/I^)) (21) 

aoGM\ 2e arctan (^^ tan (^^)) sm{^/A^ ) 



2(1 -e2) c^d J 713^ 



where 



"Gi? = ^3^ • (22) 



3GM 

i;2T 

Both the solutions Uq^jj2 and ii/^2 have the same structure of the standard solution for Schwarzschild 
backgrounds |2HI22j. the first term is a constant that can be neglected, the second term has a 
period that is a multiple of the period of solution Uq jj2 (J16p contributing a small correction to 
the orbital period and the last term monotonically grows with increasing ip contributing to the 
orbital precession. This last result is justified by noting that the analytic continuation of the 
inverse of a function corresponds to the argument of the function (in this way arctan(tanc^) = if 
increases monotonically with cp). Due to the corrections to the Keplerian orbit's solution be- 
ing small when compared to the dominant term, we can expand the trigonometric functions to 
lower order [2ll|22]: cos{VA(p) = 1 - A'^tf'^/2 + 0{(f^), <fsm{VAip) = iVAip"^/! + 0(99^) and 

avctan(y^{U^eY/JTT~e)tan(^/A^p/2]^ = Ay/{1 - e)/(l + 6)932/2 + 0{ip'^). Hence, neglecting 

the constant terms in the solutions "Ugr.h^ and Uff2 (|2ip and gathering the several terms and 
respective coefficients for these lower order expansions, we obtain the full solution u ([19]) 



Ay?GR 

2tt 



- ( 1 + e COS ( ( 1 — ^— j f] ] + Uosc.gr + W0SC.H2 



"GR 



- + acR oc [1 + 



2tt 2 --'^ ^ \ 2,d 



d^m 



+ (l_e)(l + e)f Wd 2(l-e2)GMj+^(^°" 



(23) 



Uosc.GK = — 771-^ cos(2(p) 
Qd 



{5c - 5A)aGV. 2 in \ ^ lu ^n2 -4 + e2 + 3e2 cos(299) rM u^\ 

-OSC.H2 = 1 e cos(2v.) + (go^) 4(l-e2)2GM(l + ecos^) +^(''°^ 

where Ac^GR/(2vr) is the standard precession per turn of the orbit due to General Relativity 
corrections on Schwarzschild backgrounds and IS.ip ^2 / {2,11) is the precession per turn of the orbit 
due to the ELA metric background. As for the factor Uqsc.gr it is the General Relativity oscillatory 
factor correction to the orbit solution obtained for Schwarzschild backgrounds and u^sc.H^ is the 
oscillatory factor correction due to the ELA metric background. 

To compute the observable period correction to the orbits due to the ELA metric background it 
is enough to consider the definition of the constant of motion J dr = —dip/v? ([9]) and integrate 



the infinitesimal proper time displacement dr over one turn of the orbit such that we obtain 



2n 



dip 



u^ 



2tt 



dip ■ 



1 



1 



2u. 



osc.gr 



2'Un 



;.H2 



(24) 



/O "■" -" ^0 "0 V ""0 ^0 

To directly compare the General Relativity corrections to the orbital period on the ELA met- 
ric background with the Keplerian orbital period and the General Relativity corrections on 
Schwarzschild backgrounds this integral can be factorized into the 3 components 

T = To + ATgr + ATh2 , 



1 

Jo 



2tt 



-r / dip-^ 



Vgm 



AT, 



GR 



+ 



AT, 



2tt 



dip 



^osc.GR _ 37r VGM rl^^^ e 



c2(l-e2)2 

i-2n 



(25) 



m 



{To + ATgr) {Sa - Sb) 



2 



dip ■ 



:.H2 







(5. 



Vgm 



STrVGMrl^^^e' 



2c2(l-e2)2 



3 

Svr {6a - Sc) rl^^^ e^ 



Vgm 



c2(l-e2) 



+ 



vr r 



9 
2 
l.orb 



(4 + 9e2) Hi 



{GMf. 



0{Ht 



where Tq is the classical Keplerian orbit period corresponding to the solution uq dZ]), ATgr is 
the standard General Relativity period correction on Schwarzschild backgrounds corresponding 
to solution Uosc.gr (|23p and ATjj2 is the General Relativity period correction on the ELA metric 
background corresponding to solution u^scH^ (|23]) . 



Planet 


ATgr (s/yr-') 


— - — - (arcsec/ century ^) 

27r 


Mercury 


+2.35 X 10-3 


10.35 


Venus 


+ 1.36 X 10-9 


5.30 


Earth 


+3.64 X 10-« 


3.84 


Mars 


+2.38 X 10-5 


2.54 


Jupiter 


+4.95 X 10-^ 


0.74 


Saturn 


+4.24 X lO-'^ 


0.40 


Uranus 


+ 1.21 X lO-'^ 


0.20 


Neptune 


+8.44 X 10-" 


0.13 


Pluto 


+5.15 X 10-5 


0.10 



Table 2: Standard General Relativity corrections to the orbital period ATqr Ii25\) and orbital preces- 
sion AipQ-n/2'K i23\) for each planet in the Solar System. These corrections are computed for the Solar 
Schwarzschild background with respect to the respective quantities for Keplerian orbits. 



The values for the standard General Relativity corrections to the orbital period and orbital preces- 
sion on Schwarzschild backgrounds are listed in table [2l These corrections correspond to orbital 
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motion on Schwarzschild backgrounds being well known and already accounted for in the model 
employed in the numerical analysis of the Solar System dynamics. 

In addition, although on Schwarzschild backgrounds the orbital radius is not varying over time, 
on expanding backgrounds such as the ones described by the ELA metric, it is expected that the 
radius does vary as the background expands. The analytical solutions computed so far do not 
allow to estimate such variation for the orbital radius as we have approximated the ELA metric 
background by a static background with fixed Hubble rate H = Hq ([8]). In the next section, 
from conservation of angular momentum, we estimate the orbital radius variation by considering 
approximately circular orbits. 

3 Circular Orbits Approximation: Time Varying Orbital Radius 

In this section, with the objective of estimating the orbital radius variation on backgrounds de- 
scribed by the ELA metric, we are analysing circular orbits on such backgrounds. In the non- 
relativistic velocity limit and for relatively small values of the radial coordinate (ri <^ Ih = c/H) 
the radial acceleration is 

h - -c^r% « -^ + ^-^ + F^. + Oir! H') , (26) 

rf c^ rf 

where in the right hand side of the equation, the first term is the usual classical Newton gravita- 
tional acceleration, the second term is the standard General Relativity correction on Schwarzschild 
backgrounds and the third term is the General Relativity correction of order H^ for backgrounds 
described by the ELA metric ^ 



2GMY'' f (l-ao)GM / 2GM\^-^" 



To derive the orbital velocity let us consider the constant of motion corresponding to conservation 
of angular momentum, J = —rf dip/ dr. Particularizing to circular orbits for which the orbital 
velocity is constant, Uorb = V—rirl, and considering the usual definition of angular velocity if = 
uj = Uorb/'^i we obtain the following definition for the angular momentum J^j^.^ = —7^ rf ri | , ^ 
such that for an orbit of radius ri.orb we obtain 

Ji, « GM n.o,b U + (^^) ' (1 - ^^^ "] - Fh^ -Lrb • (28) 

Here 7 = dt/dr is the relativistic factor for the ELA metric (jl]) and dotted quantities represent 
derivation with respect to the coordinate time t. Specifically, in the limit of non-relativistic velocity 
x^ <^ c, it is 

d'^ri 2 •• ^1 

d72" ^ ^ ''^ ^ ^ 2GM (Hri)'^(-^ 2GaP]^ ' ^^^^ 

(Pri \ c J y c^ n J 

For circular orbits, the main effect obtained due to the corrections on the ELA metric background 
correspond to a time varying radius. Such effect can be verified from conservation of angular 
momentum. To lowest order in time, H is expressed as H{t) ~ Hq — qq Hq t such that assuming 
a non-varying Gravitational constant G = and non-varying mass M = we are left with the 
only possibility of a time- varying orbital radius ri.orb = 0. Hence differentiating equation (|28p and 
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solving the equation Jdrc = for ri.orb we obtain, to lowest order in Hq, the time dependence of 
the orbital radius 



ri.orb ^. 2qoiHon,oAf (^ 2GM 



'^l.orb GM V c2 ri.orb 



"0 I 1 
2 ^2 



A f. (2-ao)GM\/ 2GM \^-^\ ^.. 



(30) 



As expected from cosmological expansion this expression increases with the orbital radius ri.orb and 
decreases with the mass M. Consistently at very large radius (ri ^ Ijj = c/H) the gravitational 
potential is negligible (1/ri ~ 0) such that pure cosmological expansion is asymptotically recovered 
and no stable orbits exist (ri/ri ~ 2qQH^rf > 0). 

As for the specific dependence of the orbital radius variation on the parameter ao it is positive 
for small values of qq ~ being of the same order of magnitude of the pure expansion effects, 
for growing positive values of this parameter, the radius variation decreases having a negative 
minimum value and then asymptotically vanishing in the limit ao ~^ +oo • This is actually 
expected, we note that in this limit the shift function is null, limQ,(,_>>-|-oo(l — f^sc) = 0, such that 
we exactly recover the SC metric, hence a Ricci flat space-time for which ri/ri is exactly null 
for all orbits. As for growing negative values of this parameter the radius variation increases up 
to a maximum positive value and then decreases monotonically. For large negative values of this 
parameter ao ^ the corrections with respect to Schwarzschild backgrounds become significantly 
higher with ri/ri < being unbounded from below. As an example of the typical values of 
ri,orb/'>'i.orb as a function of the parameter ao are plotted in figure [2] the values of ri.orb/ fi.orh for 
the Earth-Moon orbit and for Sun- Venus, Sun-Earth and Sun-Mars orbits. 

We further note that the estimative for the orbital radius variation just computed is a valid 
approximation for elliptical orbits of small eccentricity e ^ 1, hence a fairly good approximation 
for all planetary orbits in the Solar system except for Mercury's and Pluto's orbits for which 
e ~ 0.2 such that these estimative correspond at most to a rough approximation to the orbital 
radius variation for both these planets. 

Next, considering the General Relativity corrections on backgrounds described by the ELA metric, 
we will show that it is possible to map these corrections to the heuristic variation of the standard 
of measurement discussed in the introduction (either the variation of the AU or the variation of 
G), by matching the corrections to Kepler's third law on such backgrounds to the fitted value of 
the parameter AU . 

4 Variation of the Astronomical Unit 

As already discussed in the introduction the heuristic variation of the AU (or equivalently the 
variation of the Gravitational constant G) obtained from numerical analysis of the Solar System 
dynamics indicates that some unmodelled gravitational interaction is present. 

Here we are employing the ELA metric background to model such corrections to gravitational 
interactions. From the definition of the AU ([T]) and considering a non- varying Gravitational 
constant G = the effects that can match the heuristic variation of the AU are the decrease of 
the Sun's mass, the orbital period corrections and the orbital radius variation. The central mass 
reduction has the effect of increasing both the orbital radius and the orbital period such that 
angular momentum is conserved ^ and the equality ([T]) is maintained. Hence it is mapped to 
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(a) Earth-Moon 




2x10" 4x10" 6x10" 8x10" IxlO'- 




(b) Sun- Venus 




1x10** 2x10* 3x10^ 4x10* 




(d) Sun -Mars 



" 1- — 



2x10* 4x10^ 6x10* 8x10* IxlO' 



Figure 2: Examples of the profiles of the time variation rate of the orbital radius ri.orb/'^i.orb i30\) as a 
function of the parameter ao assuming the circular orbits approximation: 

(a) for the Earth-Moon orbit, the maximum positive variation is ri.orb/^'i.orb — 1.440 x 10^*^ s^^ cor- 
responding to ao = —3.084 x 10^" and the minimum negative variation for positive ao is ri.orb/''i.orb ~ 
-2.173 X IQ-'^'^s-^ corresponding to a„ = 2.396 x 10" ; 

(b) for the Sun-Venus orbit, the maximum positive variation is ^i. orb Z?"!. orb = 1-483 x 10^**^ s^^ corre- 
sponding to ao = —1.397 X lO'' and the minimum negative variation for positive ao is ri.orb/'"i.orb — 
-2.238 X lO"''^ s-^ corresponding to ao = 1.085 x 10^ ; 

(c) for the Sun-Earth orbit, the maximum positive variation is ri.orb/^'i.orb — 2.557 x 10^^" s^^ corre- 
sponding to ao = —3.609 x 10^ and the minimum negative variation for positive ao is ri.orb/^'i.orb = 
-3.859 X 10-'*2 s~i corresponding to ao = 2.803 x 10^ ; 

(d) for the Sun-Mars orbit, the maximum positive variation is ri.orb/^'i.orb — 9.043 x 10^'*'^ s^^ corre- 
sponding to ao = —5.498 x 10^ and the minimum negative variation for positive ao is ri.orb/^'i.orb = 
-1.365 X lO"''^ s^i corresponding to ao = 4.271 x 10^. 



the variation of the AU as a positive contribution approximately given by AU /AU ~ —Mq/Mq. 
The estimative for the Sun's mass variation due to radiation and matter emission is (see [UlTj and 
references therein) 

Me 



Mp, 



-6.7+1} x 10-1^ (yr-i), 



(31) 



hence, from this estimative we obtain a equivalent contribution for the variation of the AU of 



AAU 



lAU.M 



Me 

Mr. 



AU «+1.002l[!;^^4(cmyr-i) , 



(32) 



with AU evaluated in centimetre. As for the period correction ATff2, due to General Relativity 
corrections on backgrounds described by the ELA metric, it is mapped to a correction to the 
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variation of the AU given by the following difference 



^AU,AU,n. = — ^GMo^ J 

_ L ^^ /TW^+ATWgr \ n ^ \ (33) 



2 _\_ ATiy^uH2 

3 lyr TiAu.o 



AU (cmyr ) 



with AU evaluated in centimetre and Tiau.o and ATiau,h'^ corresponding to the period and the 
period variation over one revolution of the orbit evaluated both in the same units of time (e.g. 
second). The factor of yr^^ is explicitly written converting between the time units for Tiau.o 
and years such that 1 yr corresponds to one revolution of a Keplerian orbit with semi-major axis 
of length of lAU. We remark that one Keplerian Earth year differs from Tiau.o by 5.21s. In 
addition, a non-null orbital radius variation, is directly mapped to the variation of the AU such 
that for an orbital radius of ri^orh.iAU = 1 ^U we obtain 

AAUiAu.r = z AU [cm yr ) , (34) 

n.orh.iAu ^yr 

where AU is evaluated in centimetre and ri.orh.iAU /''"i.orh.iAU and Tiau.o are evaluated in the same 
units of time (e.g. second). 

Both the contributions ()33p and (j34p depend on the value of the metric parameter oq. To obtain a 
precise estimative for the values of these contributions it would be required a numerical analysis of 
the Solar System dynamics, however to further proceed analytically let us consider the approximate 
analytical estimative discussed in the previous sections. Specifically, the contribution (f33l) is 
evaluated by considering the period correction ([25|) computed for static analytical orbital solutions 
and, the contribution (I34p . is evaluated by considering the orbital radius variation ()30p computed 
for circular orbital solutions. These estimative are a fairly good approximation except for the 
orbital radius variation for the planet Mercury and the planet Pluto due to the relatively high 
eccentricity of their orbits. Nevertheless, as we are going to show next, the main contributions 
that are mapped to the variation of the AU are due to the decrease of the Sun's mass (I32p and 
to the General Relativity corrections to the orbital period on the ELA metric background (j33p 
such that the contribution due to the orbital radius variation (f34l) is negligible being lower than 
the remaining contributions by a factor of 10^^^ at Mercury's orbit and by a factor of 10^^*^ at 
Pluto's orbit, hence being negligible for the estimative obtained. For completeness of the analytical 
expressions we are keeping these contributions in the following derivations. 

We are proceeding to evaluate an analytical estimative for the variation of the AU by considering 
the orbits of the nine planets in the Solar System. We note that for each orbit, the corrections to 
Kepler's third law are proportional to the orbital radius ri.orb- To relate these corrections to the 
variation of the AU it is required to rescale the contributions from each planetary orbit to the AU 
length, hence to multiply each contribution by the factor j4C//ri.orb- Specifically, for a given orbit 
of semi-major axis ri oj-b) the respective mapped variation of the AU, for an observer at the Sun, 
is 

AU^^ut = (-^ + ^1:2£^^H^ + ^ 1 ^InlIlMiA AU , (35) 

\^ Mq ri.orb lyr 3 yr Torb.o ^orb.o J 

where Torb.o and ATorb are, respectively, the Keplerian orbital period and the General Relativity 
correction to the orbital period of each planetary orbit on the ELA metric background with respect 
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to Schwarzschild backgrounds. Tiau.o is the Keplerian period corresponding to the definition of 
the AU dl]) and in the last term the factor TiAU.o/Torh.o scales the temporal variation from each 
planet year to the AU year. For a given planetary orbit these corrections have the effect of a 
significant decrease on the variation of the AU for large positive values of ao ^ {AU < 0) and 
of a significant increase of the AU for large negative values of oq ^ {AU > 0). The dependence 
of the parameter AUorUt on the values of ao for the inner planets of the Solar System, Mercury, 
Venus, Earth and Mars is plotted in figure O 





(a) Mercury'-s Orbil conlribution to AU 


15000 


\ 


- 


5000 



-5000 


v_ 


- 




^ 



(b) Venus' Orbit contribution to AU 




(c) Earth's Orbit contribution to AU 



(d) Mars's Orbit contribution to AU 




Figure 3: Corrections AU orhu i35\} to the variation of the AU as a function of the constant metric param- 
eter ao for (a) the orbit of Mercury; (h) the orbit of Venus; (c) the orbit of Earth; (d) the orbit of Mars. 



As for the corrections to Kepler's third law as perceived for Earth based range measurements it 
is computed by evaluating the difference between the corrections corresponding to the geodesic 
motion of Earth and the corrections corresponding to the geodesic motion of the planet for which 
the range measurement is being considered. Hence, for range measurements between Earth and 
any other planet in the Solar System, these corrections are mapped to a variation of the AU by 



AU 



range 



'"l. orb. Earth ■^'~' orbit. Earth '"l. orb. planet ■"^'-^ orbit. planet 



1 



(A'^l. Earth-planet) 



(36) 



where AC/orbit. Earth and ^C/orbit.pianet Correspond to AC/orbit (|35]) evaluated for the orbit of Earth 
and the orbit of each planet, respectively. The distance (Ari, Earth-planet) corresponds to the aver- 
age distance between Earth and the planet assuming that both have Keplerian orbits. Planetary 
orbits of distinct bodies in the Solar System are generally not locked to each other, however to 
compute an estimative to this average value we are parameterizing Earth's orbit by the angular 
parameter Tpianet v/(l s) and the planet orbit, for which the range measurement is being consid- 
ered, by TEarth </'/(! s) with if G [0, 27r[ such that the average is performed over T = Tpianet/ {'^ s) 
periods of the Earth orbit and Te = TEarth/ {^s) periods of the planet orbit. Further aligning 
both the Earth and planet perihelion with the x axis, considering a rotation of the planet orbit 
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around the x axis by its orbital inclination angle 6orb followed by a rotation around the z axis 
by the angle between the longitudes of the planet's perihelion and Earth's perihelion 6per — ^per.E 
(given in table [1]), we obtain the following Keplerian estimative for the average distance between 
Earth and any other planet in the Solar System 



(Ari.Earth-planet) = ^ / dip^/XlP^+^Ky^+Xz^ , 



. cos(T(p) 1 / ,^ , ,^ ^ , 

Ax = -— — -— COs(Te; If) COS[dper - Oper.E) 

+ sin(TE ^) COs{5orb) Sm{5per - Sper.E, 

Ay = -— — -— ( COs(Te; ^p) Sin(5per - Sper.E) 

- sin(7E (fi) COs{5orb) COs{5per - Sper.E, 

sin(7E ip) sm{6orb) 



(37) 



Az 



uq{Tev) 



where the functions uqe and uq correspond, respectively, to the inverse of the orbital radius for 
Earth Keplerian orbit and the planet Keplerian orbit ([7]). 

Further noting that, when performing numerical analysis of the Solar System dynamics, the fitted 
variation of the AU is approximately a linear effect [1] being independent of each planet's mass, 
a estimative for the average value of AU orhit dSS]) can be obtained by a simple average of the 
contributions due to each planetary orbit 

9 
/ ^ ^tV orbit. i 

Au) = ^^^- . (38) 

/ orbit 9 

As for the average value for the contribution of AC/i-ange to the experimental data measurement 
can be estimated by an average weighted by the number of events Ni divided by the respective 
rms residuals iTj 

V AU -^ 

/ ^ ^T^"-^ range. I 

(au) = '-^^ . (39) 

\ /range « ^. 

In this expression the index i runs from 1 to 9 referring to the planets in the Solar System listed in 
table [Hand AUorUtA are the variation of the AU for each planetary orbit in the Solar System ()35p . 
The values of the weights Wi = Ni/ui are computed from table 2 and table 3 of [4J being oji = 4.65, 
UJ2 = 35692.38, W3 = 0, 0^4 = 469279.08, W5 = 74.86, ojq = 348.02, 0J^ = 62.40, ojs = 64.66 and 
wg = 38.79. 

Hence, when considering the ELA metric background, there will be two distinct corrections which 
can be mapped into the heuristic fit of the parameter AU . Let us recall that the numerical 
analysis of the Solar System dynamics has, generally, two distinct procedures p3]. First the 
ephemerides are built considering as the base model only the well established General Relativity 

15 



gravitation interactions on Schwarzschild backgrounds. Then the ephemerides are numerically 
integrated by considering a wide number of parameters which generally may include corrections 
to the gravitational interactions on Schwarzschild backgrounds, for example the PPN parameters 
7, /3 and a, as well as the variation of the AU, corresponding to the parameter AU or the variation 
of the Gravitational constant corresponding to the parameter G. It is from this second numerical 
analysis that an heuristic fit to the parameter ( AU ) is obtained ([3]). Therefore when considering 

\ / Fit 

the ELA metric background to model gravitational interactions corrections to the gravitational 
interactions on Schwarzschild backgrounds it is required to consider corrections to both of these 
procedures. In particular when mapping the modelled correction to an heuristic variation of the 
AU with respect to some fixed length AUq we obtain that the average value with respect to the 
ephemerides data is {AU) ^ , = AUq + \AU) and, when fitting the ephemerides data to 

^ \ / range 

a heliocentric model of the Solar System, we obtain {AU) ^^ = AUq + \AU) — \AU) , 

^ \ / Fit \ / orbit 

where the contribution \AU\ corresponds to the heuristic fit ([3]). Hence matching these two 

\ / Fit 

expressions we obtain the following map between the heuristic fit to the parameter AU ([3]) and 
the modelled average contributions p8|) and ([39|) 



Au) = (au) + (au) . (40) 

/ Fit \ I orbit \ / range 

Recalling that the ELA metric interpolates between local gravitational backgrounds and the ex- 
panding cosmological background we note that the previous discussion is actually consistent with 
the fact that expansion effects are proportional to the distance between the observer and the 
observed event, hence for range measurements it is required to consider the difference between 
geodesical motion of the observer (at Earth) and the geodesical motion of the remaining planets 
while for heliocentric orbital motion it is considered only the geodesic motion for each planet (the 
observer is at the Sun). 

Generally we could consider a variation of the functional parameter a across the Solar System such 
that, for each planetary orbit, this metric parameter would be given by an approximate constant 
value ao.j. However an exact fit to ephemerides would require a full numerical analysis including 
the gravitational corrections due to the ELA metric background. Nevertheless, for analytical 
analysis purposes, let us simply consider an approximately constant coefficient qq across the Solar 
System. Hence, from the map (I40p and for the value of the heuristic fit ([3]) we obtain 

AU\ = (0.75 ±0.60) cm yr^i -^ ao = l-01+i:22 • (41) 

We note that the dependence of AU on oq has an infiexion point near ao = 1 (see figure [2]), the 
relatively large uncertainty on the value of ao is mainly due to the proximity to this infiexion 
point. In table [3] are listed the values of the contributions from each planetary orbit in the Solar 
System to the estimative (|4ip . The correction to the orbital period for the value of the parameter 
ao is approximately constant for all planetary orbits 



^Th2 1%^:£!M = _o.069+^-038 x 10" 



yr lorb.O 

where ATjj2 corresponds to the period correction per revolution for each planet (j25p . Although 
this correction is enough to map the heuristic fit to the variation of the AU, we note that it is 
negligible for most of other purposes, even for archaeological fits to the variation of the Solar 
System parameters we obtain at most a variation of the Earth year by ±1.7 /i over a period of 10^ 
years, hence within the uncertainty of such estimative [7|. 
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Planet 
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Earth — planet/ 

xlO"m) 


AA'C/orb 

icmyr"^ 


r 

) 


A^'C/orb.H= 

{cm yr^^) 






Af/orbit 

{cmyr~^) 


-^^ range 

(yr-1) 


Mercury 




1.552 


5.55 


X 


10' 


-21 


-omitlfA 




+1.093+^:^2} 


0.580t}:n7 


Venus 




1.700 


3.62 


X 


10- 


-20 


-ometlZ 




+1.038+1:}°} 


0.241+0:^,^3 


Earth 




- 


9.57 


X 


10- 


-20 


-0.022t}:r8 X 


10- 


-1 


+1.02411:^3 


- 


Mars 




2.525 


3.38 


X 


10- 


-19 


-O.OUtlf.t X 


10- 


-1 


+lMAt\%\\ 


-0.3G9tr2}^5 


Jupiter 




7.842 


1.35 


X 


10- 


-17 


-O.G18tl:607 X 


10- 


-1 


+i.oo4l0:}«} 


-0.8Gl+°:fo} 


Saturn 




14.28 


8.30 


X 


10- 


-17 


-0.074t^:f78 X 


10- 


-2 


+1.003iro^^ 


-OMbtlf^ 


Uranus 




28.72 


6.76 


X 


10- 


-16 


-0.026^1^2^^ X 


10- 


-2 


+1.005iro'2^ 


-0.949t°:^^? 


Neptune 




45.00 


2.60 


X 


10- 


-15 


-o.oi3tl:}^^ X 


10- 


-2 


+1.0021°:°}^ 


-0.968t°;Sf2 


Pluto 




57.23 


5.89 


X 


10- 


-15 


-0.089t^:«690 X 


10- 


-3 


+i.oo2t°:°°« 


-l.GG8+°:°f2 



Table 3: Corrections to AU mapped from the corrections to Kepler's third law on the EL A metric back- 
ground with respect to Schwarzschild backgrounds. For each planet it is listed, in the first column the average 
distance between the planet and Earth for Keplerian orbits {Ari, Earth-planet) [^^ , in the second column 
the contribution due to the orbital radius variation (the second term in equation Ii35\) ) for which the uncer- 
tainty is at least 8 orders of magnitude below the quoted values, in the third column the contribution due to 
the period correction on the EL A metric background (the third term in equation h3^(l ). in the fourth column 
the total correction to the variation of the A U for heliocentric orbital motion AU orbit i35\) and in the fifth 
column the total correction to the variation of the AU for Earth based range measurements AU range liS6]). 



As for the values for the corrections to the orbital precession and orbital radius variation for each 
planet are listed in table HI In particular the orbital precession corrections are lower by more 
than 10 orders of magnitude when compared to the respective values for Schwarzschild back- 
grounds and the orbital radius variation range from riorb/^i.orb ~ 10~^^ century~^ for mercury 
up to ri.orb/'^i.orb ~ 10~^^ century~^ , hence being well below any other estimative for these varia- 
tions [7] . We recall that a variation of the AU does not imply an orbital radius variation [3H6] and 
further remark that the orbital radial variation obtained from numerical analysis of ephemerides 
in [1] are due to the variation of the Gravitational constant G ^ 0, hence not directly comparable 
with the values obtained here, for which the quoted orbital radius variation in table H] is due to the 
ELA metric background alone. In the framework discussed here there will be an extra contribution 
to the orbital radius variation due to the decrease of the Sun's mass of order 10~^^ century"^ to 
10~^^ century~^ similarly to the fit obtained in [4[|26j. 

With respect to the ELA metric background we note that the value of the metric parameter ao 
parameterizes the local anisotropic corrections with respect to the isotropic cosmological back- 
ground, specifically space-time is locally isotropic for ao = which corresponds to the isotropic 
background described by the McVittie metric [l5j. The value of ao ~ 1-01 (j4ip corresponds to a 
relatively small perturbation to the isotropic background which, as has been shown, corresponds 
to relatively small corrections to the orbital parameters. Consistently with this discussion we re- 
mark that, when compared with a isotropic variation of the AU (or equivalently, to the variation 
of the Gravitational constant G), the corrections due to the ELA metric background to the or- 
bital period are relatively more relevant than the corrections to the orbital precession and orbital 
radius variation. It is due to the background anisotropy between the radial direction and angular 
directions that such effect is attainable. In addition we recall that, when considering point-mass 
objects, as we approach the SC horizon the metric exponent a should be greater or equal to 
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Planet 


1^" {arcsec 


century i) 


''^-°"'' (century i) 


Mercury 


-2.28 X 


10-19 


3.71 X 


10-32 


Venus 


-5.82 X 


10-19 


2.42 X 


10-31 


Earth 


-9.46 X 


10-19 


6.40 X 


10-31 


Mars 


-1.78 X 


10-18 


2.26 X 


10-30 


Jupiter 


-1.12 X 


10-17 


9.01 X 


10-29 


Saturn 


-2.79 X 


10-17 


5.55 X 


10-28 


Uranus 


-7.96 X 


10-17 


4.52 X 


10-27 


Neptune 


-1.56 X 


10-16 


1.74 X 


10-26 


Pluto 


-2.35 X 


10-16 


3.94 X 


10-26 



Table 4: Corrections to precession and orbital radius variation due to the General Relativity corrections 
for the background described by the ELA metric. The estimative uncertainty for each of these values is 
several orders of magnitude below the quoted values. 



o(^i.sc) = 3 to ensure that space-time is singularity free at this horizon. This requirement is not 
absolutely necessary as the real Sun is not a point-like mass being instead an extended spheroid, 
hence without an event horizon. Nevertheless we further note that the uncertainty on the esti- 
mative of the constant ao is relatively large and we may expect that a functional parameter a 
varying across the Solar System would possible allow for a better fit to ephemerides, hence we 
may conjecture that its value should be decreasing with growing heliocentric distances being close 
to a = 3 near the Sun. This discussion is not conclusive being required a numerical analysis of the 
Solar System dynamics including the corrections due to the ELA metric background to actually 
verify if such a profile for the values of a is the best fit to planetary motion. 



5 Conclusions 



In this work we have mapped the General Relativity corrections to Kepler's third law on back- 
grounds described by the ELA metric (j4]) to the heuristic variation of the AU estimated from 
numerical analysis of the Solar System dynamics on Schwarzschild backgrounds. These correc- 
tions together with the decrease of the Sun's mass by radiation emission fully account for the fitted 
value for the variation of the AU ([3]) . Reflecting the anisotropic nature of the ELA metric back- 
ground, the more relevant contribution to such modelling is due to the orbital period corrections, 
being the contributions from the orbital radius variation negligible. 

The constant value for the metric parameter that matches the quoted variation of the AU is 
ao = +1.01]1^^'22 (BH)) hence relatively close to the value oq = which corresponds to the isotropic 
background described by the McVittie metric. For completeness let us further note that other 
effects which may be relevant on backgrounds described by the ELA metric such as the corrections 
to the Doppler shift for range measurements are, for this range of values of the inetric parameter 
negligible, being of the same order of magnitude of the effects attributed to isotropic cosmological 
expansion |20tl27j. Following the same arguments we further conclude that, within the Solar 
System the contribution to the cosmological mass-energy density within the Solar System due to 
the ELA metric background is negligible [28] (for further details see |20]). 
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Hence we have shown that the heuristic variation of the AU (or equivalently the heuristic vari- 
ation of the Gravitational constant G) OH] can alternatively be modelled by the ELA metric 
background parameterizing the corrections to gravitational interactions within the Solar System. 
For analytical analysis purposes we have considered a constant metric parameter oq for which we 
obtain a relatively high uncertainty. More generally a radially symmetric functional parameter 
with varying value across the Solar System would, in principle, allow for a significant reduction 
of such uncertainty. The results obtained here are enough to motivate a numerical analysis of 
the Solar System dynamics including the ELA metric background corrections to gravitational 
interactions. This framework solves the problem of the unwelcome variation of the measurement 
standard (whether the AU whether G) and constitutes a playground for testing the ELA metric 
background in the most well known of all the astrophysical systems, the Solar System. We leave 
such study to another work. 
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